Sterile neutrinos and Rk 



A. Vicente 

Laboratoire de Physique Theorique, CNRS - UMR 8627, Universite de Paris-Sud 11, F-91405 
Orsay Cedex, France 

E-mail: avelino .vicente@th.u-psud.fr 

Abstract. We study the violation of lepton flavour universality in light meson decays due to 
the presence of non-zero mixings between the active neutrinos with new sterile states. The 
modified Wiu vertices, arising from a non-unitarity leptonic mixing matrix intervening in 
charged currents, might lead to a tree-level enhancement of Rp = r(P — >■ eiy)/r{P — >■ fiv), 
with P — K,TT. These enhancements are illustrated in the case of the inverse seesaw, showing 
that one can saturate the current experimental bounds on Arx (and Ar^), while in agreement 
with the different experimental and observational constraints. 



1. Introduction 

Lepton flavour universality (LFU) is a distinctive feature of the Standard Model (SM). The 
different lepton families couple with exactly the same strength to the gauge bosons. This leads 
to concrete predictions in electroweak precision tests, which can only distinguish among lepton 
families by the different charged lepton masses. Any deviation from the expected SM theoretical 
results would signal the presence of New Physics (NP). In this work we concentrate on light 
meson {K and tt) leptonic decays which, in view of the expected experimental precision, have a 
unique potential to probe deviations from the SM regarding lepton universality. 

In the SM, the dominant contribution to T(P — )• ii^) (P = K,Tr) arises from W boson 
exchange. One may be afraid about potential hadronic uncertainties; however, by considering 
the ratios 

^ T{K+ ^ e+z.) ^ r(7r+ ^ e+t^) 

^ " T{K+ ^ n+u) ' " " r(^+ ^ fi+u) ' ^ ' 

the hadronic uncertainties are expected to cancel out to a good approximation. In order to 
compare the experimental bounds with the SM predictions, it is convenient to introduce a 
quantity, Arp, which parametrizes deviations from the SM expectations: 

i?^^P = i?SM ^ equivalently Arp = - 1 . (2) 

The comparison of theoretical analyses [li [2j with the recent measurements from the NA62 
collaboration [3j and with the existing measurements on pion leptonic decays [1] 

Rj^ = (2.477 ± 0.001) X lO'^ , = (2.488 ± 0.010) x 10"^ , (3) 

^SM ^ (1.2354 ± 0.0002) x 10-^ i?^^P = (1.230 ± 0.004) x 10"^ (4) 



suggests that observation agrees at la level with the SM's predictions for 

ArK = (4 ±4) X 10"^ Ar^ = (-4 ±3) x 10"^ (5) 

The current experimental uncertainty in Ar^ (of around 0.4%) will be further reduced in the 
near future, as one expects to have 6Rk/Rk ~ 0.1% [5j, which can translate into measuring 
deviations Ar^ ~ 0(10"'^). Similarly, there are also plans for a more precise determination of 
Ar, 16, 7J. 

New contributions to Arp have been extensively discussed in the literature, especially in the 
framework of models with an enlarged Higgs sector. In the presence of charged scalar Higgs, new 
tree-level contributions are expected. However, as in the case of most Two Higgs Doublet Models 
(2HDM), or supersymmetric (SUSY) extensions of the SM, these new tree-level corrections 
are lepton universal [8j. In SUSY models, higher order non-holomorphic couplings can indeed 
provide new contributions to Rp O UHl IHl III [l3] , but in view of current experimental bounds 
(collider, i?-physics and r-lepton decays), one can have at most Arx < 10~^ in the framework 
of unconstrained minimal SUSY models Corrections to the Wiu vertex can also induce 

violation of LFU in charged currents. However, if these appear at the loop level, as referred to 
in [TU], the effect is expected to be of order (a/47r) x (m^/A^p) (Anp being the new physics 
scale), generally well below experimental sensitivity. 

Here we consider a different alternative: The tree-level corrections to charged current 
interactions once neutrino oscillations are incorporated into the SM \15\ I16j. In this case, and 
working in the basis where the charged lepton mass matrix is diagonal, the flavour-conserving 
term oc glj^^^PLi'jW~ now reads 

- = -^Ul%rPL^^W;: + C.C. , (6) 

where Ul^ is a generic leptonic mixing matrix, i = 1, . . . ^n^ denoting the physical neutrino states 
(not necessarily corresponding to the three left-handed SM states = vl) and j = 1, ... ,3 the 
charged lepton flavour. In the case of three neutrino generations, Ul^ corresponds to the unitary 
PMNS matrix and flavour universality is preserved in meson decays: since one cannot tag the 
flavour of the final state neutrino (missing energy), the meson decay amplitude is proportional 
to {UyUl)jj = 1, and thus no new contribution to Rp is expected. However, in the presence of 
sterile states, the W^l' vertex is proportional to a rectangular 3 x matrix Ul^ , and the mixing 
between the left-handed leptons vlAl corresponds to a 3 x 3 block of Ul^ , 

C^PMNS — ^ ^^PMNS = (1 — f?) f^PMNS • (7) 

The larger the mixing between the active (left-handed) neutrinos and the new states, the more 
pronounced the deviations from unitarity of f/pMNS) parametrized by the matrix r/ [14J. The 
active-sterile mixings and the departure from unitarity of f/pMNS can be at the source of the 
violation of LFU in different neutrino mass models which introduce sterile fermionic states to 
generate non-zero masses and mixings for the light neutrinos [16] . 

Corrections to the Wlv vertex can arise in several scenarios with additional (light) singlet 
states, as is the case of z/SM [1^, low-scale type-I seesaw [18] and the Inverse Seesaw (ISS) [19], 
among other possibilities. This clearly shows the potentiality of the mechanism under discussion, 
which can be present in many different models. 

In the next section we provide a model- independent computation of Arp in the presence 
of additional fermionic sterile states; we then briefly review in Section [3] the most important 
experimental and observational constraints on the mass of the additional singlet states. In 
Section we consider the case of the inverse seesaw to give a numerical example of the impact 
of sterile neutrinos on Arp. Our concluding remarks are summarised in Section [5j 



2. Arp in the presence of sterile neutrinos 

Let us consider the SM extended by Ng additional sterile states. The matrix element for the 
meson decay P — )• IjVi can be generically written as 

M^, = u,^ {A'^Pr + B'^Pl)vi^ . (8) 

No sum is implied over the indices of the outgoing leptons Notice that i = 1, . . . ,3+Ns. The 
expressions for A and B can be easily obtained from the usual 4-fermion effective hamiltonian 
obtained after integrating out the W boson in Eq. These are 

(Af = (A'^f = -4 Gf FcIm fp Ul' * mi^ ; (9) 
{Bf^ = {B'^r = 4 Gf ^ckm fp Ui' * m,, , (10) 

where fp denotes the meson decay constant and rai.^y. the mass of the outgoing leptons. 
The expression for Rp is finally given by 

RP = ^^^^^ ^ith (11) 



F'^ = \W;\^ and G'^ 



2/2, 2 \ / 2 2\2 



f 2 2 2 \2 ^ 2 2 

[mp - nii^ -m^J - 4m;^.m^^ 



1/2 

(12) 



The result of Eq. (jlip has a straightforward interpretation: F^^ represents the impact of new 
interactions (absent in the SM), whereas G*-' encodes the mass-dependent factors. The SM result 
can be easily recovered from Eq. ()lip . in the limit rrii,- = and Ul^ = Sji, 

J.SM _ ml (ml - mlf 

-ftp — 9/2 9^9 ' V 

[rrip — m^j^ 

to which small electromagnetic corrections should be added pp. 

Using the results in Eqs. ([TT]) and (fT3|) . we obtain a general expression for Arp 

mliml, - mlY y^^^ix pmi Qmi 

— A''- P M-' L^m=l ^ (14) 

ml{m], - 771.2)2 ^Wi|x pm2 Qn2 

Thus, depending on the masses of the new states (and their hierarchy) and most importantly, 
on their mixings to the active neutrinos, Arp can considerably deviate from zero. In order to 
illustrate this, we consider two regimes: 

• Regime (A): All sterile neutrinos are lighter than the decaying meson, but heavier than 
the active neutrino states, i.e. m^^^"^ <^ rriy^ ^ ^P 

• Regime (B): All sterile neutrinos are heavier than mp 

Notice that in case (A) , all the mass eigenstates can be kinematically available and one should 
sum over all S + Ng states; furthermore there is an enhancement to Arp arising from phase space 
factors, see Eq. ([12]). 



3. Constraints on sterile neutrinos 

We review in this section the experimental and observational bounds on the mass regimes and 
on the size of the active-sterile mixings that must be satisfied. 

First, it is clear that present data on neutrino masses and mixings [20j should be accounted 
for. Second, there are robust laboratory bounds from direct sterile neutrinos searches [2H I22j. 
since the latter can be produced in meson decays such as vr^ — t- ij.^u, with rates dependent on 
their mixing with the active neutrinos. Negative searches for monochromatic lines in the muon 
spectrum can be translated into bounds for tti^^ — combinations, where parametrizes 
the active-sterile mixing. The non-unitarity of the leptonic mixing matrix is also subject to 
constraints. Bounds on the non-unitarity parameter rj (Eq. dZD), were derived using Non- 
standard Interactions [23j; although not relevant in case (A), these bounds will be taken into 
account when evaluating scenario (B). 

The modified Wiu vertex also contributes to lepton flavour violation (LFV) processes. The 
radiative decay // — )• ey, searched for by the MEG experiment [24J, is typically the most 
constraining observable \j. The rate induced by sterile neutrinos must satisfy [311 [32] 

^^(^ ^ = ^^4^\H,e\' < 2.4 X 10-12 , (15) 

2 

where H^e = Y.iUfUl"G^{^^), with the loop function and the mixing matrix 

defined in Eq. ([Hj). Similarly, any change in the W^v vertex will also affect other leptonic 
meson decays, in particular B ^ iv; the following bounds were enforced in the analysis: 
BR{B ev) < 9.8 x 10"^ BR{B ^v) < 10"^ and BR{B tv) = (1.65 ± 0.34) x lO""^ [33|. 

Important constraints can also be derived from LHC Higgs searches [?] and electroweak 
precision data |35j . They will also be considered in our numerical analysis. 

Under the assumption of a standard cosmology, the most constraining bounds on sterile 
neutrinos stem from a wide variety of cosmological observations [361 122] . These include Large 
Scale Structure data, X-ray searches (which can be produced in z/j — >■ z^j7), Lyman-a limits, the 
existence of additional degrees of freedom at the epoch of Big Bang Nucleosynthesis and Cosmic 
Microwave Background data. However, all the above cosmological bounds can be evaded if a 
non-standard cosmology is considered. In fact, the authors of Ref. [37J showed that the above 
cosmological constraints disappear in scenarios with low reheating temperature. Therefore, we 
will allow for the violation of the latter bounds, explicitly stating it. 

4. A numerical example: Ar/< in the inverse seesaw 

Although the generic idea explored in this work applies to any model where the active neutrinos 
have sizeable mixings with some additional singlet states, we consider the case of the Inverse 
Seesaw [TD] to illustrate the potential of a model with sterile neutrinos regarding tree-level 
contributions to light meson decays. As mentioned before, there are other possibilities [TTl [T^ . 

^.1. The inverse seesaw 

In the ISS, the SM particle content is extended by generations of right-handed (RH) neutrinos 
vji and nx generations of singlet fermions X with lepton number L = —1 and L = -|-1, 
respectively p^£j (such that nn + nx = Ng). In our numerical application we will focus on 

^ Recently, it has been also noticed that in the framework of low-scale seesaw models, the expected future 
sensitivity of p — e conversion experiments can also play a relevant role in detecting or constraining sterile 
neutrino scenarios [251 1261 [?7ll28j . This is also the case in the supersymmetric version of these models, even when 
the sterile neutrinos are heavier [291 130j . 



the case ur = nx = 3. The lagrangian is given by 



Ciss = CsM + Y^^URiLjH + Mr,^ umXj + ]^^xxijX'iXj + h.c. (16) 

where i,j = 1,2,3 are generation indices and H = ia2H* . Notice that the present lepton 
number assignment, together with L = +1 for the SM lepton doublet, implies that the "Dirac"- 
type right-handed neutrino mass term Mr., conserves lepton number, while the "Major ana" 
mass term ^Xij violates it by two units. 

The left-handed neutrinos mix with the right-handed ones after electroweak symmetry 
breaking. This leads to an effective Majorana mass for the active (light) neutrinos. Assuming 
fix *C rriD ^ Mr, where rriD = "^^i'^' with v the vacuum expectation value of the SM Higgs 
boson, one obtains 

rriu ^ m^MJi fixMj^^mn . (17) 

The remaining 6 sterile states have masses approximately given by M^, ~ Mr. Small corrections 
can be added to these results, but they are typically negligible [38] . 

In what follows, and without loss of generality, we work in a basis where Mr is a diagonal 
matrix (as are the charged lepton Yukawa couplings). Yi, can be written using a modified 
Casas- Ibarra parametrisation [39j (thus automatically complying with light neutrino data), 

/2 

Yu = ^vWMRy^ul^^s^ (18) 

where \/m^ is a diagonal matrix containing the square roots of the three eigenvalues of rrii, 
(cf. Eq. ()17p ): likewise \/m is a (diagonal) matrix with the square roots of the eigenvalues of 
M = Mrh'^M]^. V diagonahzes M as YMV"^ = M, and R is a 3 X 3 complex orthogonal 
matrix, parametrized by 3 complex angles, encoding the remaining degrees of freedom. 

The nine neutrino mass eigenstates enter the leptonic charged current through their left- 
handed component (see Eq. ([6]), with i = 1, . . . , 9, j = 1, . . . , 3). The unitary leptonic mixing 
matrix C/^ is now defined as U^MUy = diag(mj). Notice however that only the rectangular 
3x9 sub-matrix (first three columns of Uu) appears in Eq. ([6]) due to the gauge-singlet nature 
of ur and X. 

4-2. Numerical evaluation of Ar^ in the inverse seesaw 

We numerically evaluate the contributions to Rk in the framework of the ISS and address the two 
scenarios discussed before, which can be translated in terms of ranges for the (random) entries 
of the Mr matrix: regime (A) {m^^ < mp) - Mr^ G [0.1,200] MeV; regime (B) {m^^ > mp) - 
Mri G [1, 10^] GeV. The entries of nx have also been randomly varied in the [0.01 eV, 1 MeV] 
range for both cases. 

The adapted Casas- Ibarra parametrisation for Yy, Eq. (jl8p . ensures that neutrino oscillation 
data is satisfied (we use the best-fit values of the global analysis of Ref. [20j and set the 
CP violating phases of C/pMNS to zero). The R matrix angles are taken to be real (thus no 
contributions to lepton electric dipole moments are expected), and randomly varied in the range 
Oi E [0,27r]. We have verified that similar Ar^^ contributions are found when considering the 
more general complex R matrix case. 

In Figs. [1] we collect our results for Arx in scenarios (A) - left panel - and (B) - right 
panel, as a function of fj, which parametrizes the departure from unitarity of the active neutrino 
mixing sub-matrix f/pMNSi = 1 — |Det(C/pMNs)|- Although the cosmological constraints are 
not always satisfied, we stress that all points displayed comply with the different experimental 
and laboratory bounds discussed before. For the case of scenario (A), one can have very large 
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Figure 1. Contributions to Arx in the inverse seesaw as a function of fy = 1 — |Det(C/pMNs)|' 
regimes A (left) and B (right). The upper (lower) dashed line denotes the current experimental 
limit (expected sensitivity). On the right panel, red points denote cases where > W'^. All 
points comply with experimental and laboratory constraints. Points in (B) are also in agreement 
with cosmological bounds, while those in (A) require considering a non-standard cosmology. 



contributions to Rk, which can even reach values Arx ~ C(l) (in some extreme cases we 
find Arx as large as ~ 100). The hierarchy of the sterile neutrino spectrum in case (A) is 
such that one can indeed have a significant amount of LFU violation, while still avoiding non- 
unitarity bounds. Although this scenario would in principle allow to produce sterile neutrinos 
in light meson decays, the smallness of the associated Y^, (< 0(10"^^)), together with the loop 
function suppression (G-y), precludes the observation of LEV processes, even those with very 
good associated experimental sensitivity, as is the case of /x — )• ej. The strong constraints from 
CMB and X-rays would exclude scenario (A); in order to render it viable, one would require a 
non-standard cosmology. 

Despite the fact that in case (B) the hierarchy of the sterile states is such that non- 
unitarity bounds become very stringent (since the sterile neutrinos are not kinematically viable 
meson decay final states), sizeable LFU violation is also possible, with deviations from the SM 
predictions again as large as Ar^ ~ C(l)- Although one cannot produce sterile states in meson 
decays in this case, the large Y^, open the possibility of having larger contributions to LFV 
observables so that, for example, BR(/i — )• 67) can be within MEG reach. 

Although we do not explicitly display it here, the prospects for Ar^^ are similar: in the same 
framework, one could have Ar^^ ~ 0{ArK), and thus Arjr ~ in both scenarios. Depending 
on the singlet spectrum, these observables can also be strongly correlated: if all the sterile states 
are either lighter than the pion (as it is the case of scenario (A)) or then heavier than the kaon, 
one finds Ar^ ~ Ar^- This is a distinctive feature of our mechanism. 

5. Concluding remarks 

The existence of sterile neutrinos can potentially lead to a significant violation of lepton flavour 
universality at tree-level in light meson decays. As shown in this study, provided that the active- 
sterile mixings are sufficiently large, the modified Wiv interaction can lead to large contributions 
to lepton flavour universality observables, with measurable deviations from the standard model 
expectations, well within experimental sensitivity. This mechanism might take place in many 
different frameworks, the exact contributions for a given observable being model-dependent. 

As an illustrative (numerical) example, we have evaluated the contributions to Rk in the 
inverse seesaw extension of the SM - a truly minimal extension of the SM - , for distinct 



hierarchies of the sterile states. Our analysis reveals that very large deviations from the 
SM predictions can be found (Arj^ ~ C'(l)) - or even larger, well within reach of the NA62 
experiment at CERN. This is in clear contrast with other models of new physics (for example 
unconstrained SUSY models, where one typically has Arx < 0(10-^)). 
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